MACAULAY 2 AND THE GEOMETRY OF SCHEMES 



GREGORY G. SMITH AND BERND STURMFELS 

This tutorial illustrates how to use Dan Grayson and Mike Still- 
man's computer algebra system, Macaulay2, to study schemes. The 
examples are taken from the homework for an algebraic geometry class 
given at the University of California, Berkeley in the fall of 1999. This 
graduate course was taught by the second author with assistance from 
the first author. Our choice of problems follows the material in David 
Eisenbud and Joe Harris' textbook The Geometry of Schemes. In fact, 
four of our ten problems are taken directly from their exercises. 

Distinguished open sets. We begin with a simple example involv- 
ing the definition of an affine scheme; see section 1.1.4 in Eisenbud 
and Harris (1999). This example also indicates some of the subtleties 
involved in working with arithmetic schemes. 

Question 1. Let R = 7*[x,y, z] and X = Spec(-R); in other words, X 
is affine 3-space over the integers. Let f — x and consider the basic 
open subset Xf C X . 

(a) If ei = x+y+z, e 2 = xy+xz+yz ande 3 = xyz are the elementary 
symmetric functions then the set {X e .} 1<i<3 is an open cover of 
Xj. 

(b) If pi = x + y + z, p 2 = x 2 + y 2 + z 2 and p 3 = x 3 + y 3 + z 3 are the 
power sum symmetric functions then {X p -} 1<i<3 is NOT an open 
cover of Xf. 

Solution. By Lemma 1-6 in Eisenbud and Harris (1999), it suffices to 
show that e±, e 2 and e% generate the unit ideal in Rf. This is equivalent 
to showing that x m belongs to the i?-ideal (ei, e 2 , e%) for some m G N. 
In particular, the saturation ((ei,e2,e3) : x°°) is the unit ideal if and 
only if {X ei } 1<i<3 is an open cover of Xf. Macaulay2 allows us to 
work with homogenous ideals over Z and we obtain: 

11 : R = ZZ[x, y, z] ; 

12 : elementaryBasis = ideal (x+y+z, x*y+x*z+y*z, x*y*z) ; 

02 : Ideal of R 

13 : saturate (elementaryBasis , x) 

03 = ideal 1 

o3 : Ideal of R 
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Similarly, to prove that {X Pi } 1<i<3 is not an open cover of Xf, it 

is enough to show that ((pi,P2,Ps) '■ is not the unit ideal. We 
compute this saturation: 

14 : powerSumBasis = ideal (x+y+z, x~2+y~2+z~2, x~3+y~3+z~3) ; 

04 : Ideal of R 

15 : saturate (powerSumBasis , x) 

2 2 

05 = ideal (6, x + y + z, 2y - y*z + 2z , 3y*z) 

o5 : Ideal of R 

However, working over the field Q, we find that ((pi,P2,P3) '■ x °°) is 
the unit ideal. 

16 : S = QQ[x, y, z] ; 

17 : powerSumBasis = ideal(x+y+z, x~2+y~2+z~2, x~3+y~3+z~3) ; 

07 : Ideal of S 

18 : saturate (powerSumBasis , x) 

08 = ideal 1 

08 : Ideal of S 

Irreducibility. The study of complex semisimple Lie algebras gives 
rise to an important family of algebraic varieties called nilpotent or- 
bits. To illustrate one of the properties appearing in section 1.2.1 of 
Eisenbud and Harris (1999), we examine the irreducibility of a partic- 
ular nilpotent orbit. 

Question 2. Let X be the set of complex 3x3 matrices which are 
nilpotent. Show that X is an irreducible algebraic variety. 

Solution. A 3 x 3 matrix M is nilpotent if and only if its minimal 
polynomial divides T fe , for some k e N. Since each irreducible factor 
of the characteristic polynomial of M is also a factor of the minimal 
polynomial, we conclude that the characteristic polynomial of M is T 3 . 
It follows that the coefficients of the characteristic polynomial (except 
for the leading coefficient which is 1) of a generic 3x3 matrix define 
the algebraic variety X. 

To show X is irreducible over C, it is enough to construct a rational 
parameterization of X; see Proposition 4.5.6 in Cox, Little, and O'Shea 
(1996). To achieve this, observe that GL n (C) acts on X by conjugation. 
Jordan's canonical form theorem implies that there are exactly three 
orbits; one for each of the following matrices: 



N = 







,N 1 = 


10 



and = 


10 
1 














Each orbit is defined by a rational parameterization, so it suffices to 
show that the closure of the orbit containing N 2 is the entire variety 
X. In Macaulay2, this calculation can be done as follows: 
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11 : S = QQ[t, y_0..y_8, a..i, MonomialOrder => Eliminate 10]; 

12 : N2 = (matrix {{0 , 1 ,0} , {0 ,0 , 1} , {0 ,0 ,0}}) ** S 

o2 = {0} I 1 I 
{0} I 1 I 
{0} I I 

3 3 

02 : Matrix S < S 

13 : G = genericMatrix(S, y_0, 3, 3); 

3 3 

03 : Matrix S < S 

To determine the entries in GN 2 G~ 1 , we use the classical adjoint to 
construct the inverse of the generic matrix G. 

14 : adj = (G,i,j) -> ( 

n := degree target G; 
m := degree source G; 

(-l)~(i+j)*det(submatrix(G, {0..(i-l), (i+1) . . (n-l)>, 
{0..(j-D, (j+l)..(m-l)})) 

); 

15 : classicalAdjoint = (G) -> ( 

n := degree target G; 

matrix table(n, n, (i, j) -> adj(G, j, i)) 

); 

16 : numerators = G*N2*classicalAdjoint(G) ; 

3 3 
06 : Matrix S < S 

17 : D = det(G) ; 

18 : M = genericMatrix(S, a, 3, 3); 

3 3 

08 : Matrix S < S 

The entries in GN 2 G~ Y give a rational parameterization of the orbit 
generated by N 2 . We give an "implicit representation" of this variety 
by using elimination theory; see section 3.3 in Cox, Little, and O'Shea 
(1996). 

19 : elimldeal = minors(l, (D*id_(S~3) ) *M-numerators) + ideal(l - D*t) ; 

09 : Ideal of S 

110 : closureOf Orbit = ideal(selectInSubring(l , gens gb elimldeal)); 

010 : Ideal of S 

Finally, we check that the closure of this orbit is equal to X scheme- 
theoretically. 

111 : X = ideal submatrix( (coefficients ({0} , 

det(M - t*id_(S~3))))_l, {1,2,3}) 

011 = ideal (a + e + i, b*d - a*e + c*g + f*h - a*i - e*i, 

- c*e*g + b*f*g + c*d*h - a*f*h - b*d*i + a*e*i) 

011 : Ideal of S 

112 : closureOfOrbit == X 

012 = true 
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Singular Points. Section 1.2.2 in Eisenbud and Harris (1999) provides 
the definition of a singular point of a scheme. In our third question, 
we study the singular locus of a family of elliptic curves. Section V.3 
in Eisenbud and Harris (1999) also contains related material. 

Question 3. Consider a general form of degree 3 in Q[x, y, z\; 

F = ax 3 + bx 2 y + cx 2 z + dxy 2 + exyz + fxz 2 + gy 3 + hy 2 z + iyz 2 + jz 3 . 

Give necessary and sufficient conditions in terms of a, . . . ,j for the 
cubic curve Proj (Q[x, y, z]/(F)) to have a singular point. 

Solution. A time consuming elimination gives the degree 12 polynomial 
which defines the singular locus of a general form of degree 3. This can 
be done in Macaulay2 as follows. We have not displayed the output 
06, as this discriminant has 2040 terms in the 10 variables a, . . . ,j. 

11 : S = QQ[x, y, z, a..j, MonomialOrder => Eliminate 2]; 

12 : F = a*x~3+b*x~2*y+c*x"2*z+d*x*y"2+e*x*y*z+ 

f*x*z~2+g*y~3+h*y~2*z+i*y*z~2+j*z~3; 

13 : partials = submatrix(jacobian matrix{{F}} , {0..2}, {0}) 

o3 = {1} I 3x2a+2xyb+y2d+2xzc+yze+z2f I 
{1} I x2b+2xyd+3y2g+xze+2yzh+z2i I 
{1} I x2c+xye+y2h+2xzf+2yzi+3z2j I 

3 1 

03 : Matrix S < S 

14 : singularities = ideal (partials) + ideal(F); 

04 : Ideal of S 

15 : elimDiscr = ideal selectlnSubringd , gens gb singularities); 

05 : Ideal of S 

16 : elimDiscr = substitute (elimDiscr, {z => 1}) ; 

06 : Ideal of S 

There is also a simple and more useful determinantal formula for this 
discriminant. It is a specialization of the formula (2.8) in section 3.2 
in Cox, Little, and O'Shea (1998): 

17 : hessian = det submatrix(jacobian ideal partials, {0..2}, {0..2}); 

18 : A = (coef f icients ({0 , 1 ,2} , submatrix( 

jacobian matrix{{F}} , {0..2>, {0})))_1; 

3 6 

08 : Matrix S < S 

19 : B = (coefficients({0,l,2}, submatrix( 

jacobian matrix{{hessian}} , {0..2}, {0>)))_1; 

3 6 

09 : Matrix S < S 

110 : detDiscr = ideal det (A I I B) ; 

010 : Ideal of S 

111 : detDiscr == elimDiscr 

011 = true 
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Fields of Definition. Schemes over non-algebraically closed fields 
arise in number theory. Our solution to Exercise II-6 in Eisenbud and 
Harris (1999) indicates one technique for working over a number field 
in MACAULAY2. 

Question 4. An inclusion of fields K L induces a map A r [ — > A^-. 
Find the images in Aq of the following points of A^ under this map. 

(a) ( X -V2,y-V2); 

(b) (x-V2,y-V3); 

(c) (x — (, y — C" 1 ) where ( is a 5-th root of unity ; 

(d) (V2x-VSy); 

(e) {V2x-V3y-1). 

Solution. The images can be determined by (1) replacing coefficients 
not belonging to K with indeterminates, (2) adding the minimal poly- 
nomials of these coefficients to the given ideal in A^ and (3) eliminating 
the new indeterminates. Here are the five examples: 

11 : S = QQ[a, b, x, y, MonomialOrder => Eliminate 2]; 

12 : la = ideal(x-a, y-a, a~2-2) ; 

02 : Ideal of S 

13 : ideal selectlnSubringd, gens gb la) 

2 

03 = ideal (x - y, y - 2) 

03 : Ideal of S 

14 : lb = ideal(x-a, y-b, a~2-2, b~2-3) ; 

04 : Ideal of S 

15 : ideal selectlnSubringd , gens gb lb) 

2 2 

05 = ideal (y - 3, x - 2) 

05 : Ideal of S 

16 : Ic = ideal(x-a, y-a~4, a~4+a~3+a~2+a+l) ; 

06 : Ideal of S 

17 : ideal selectlnSubringd, gens gb Ic) 

2 2 3 2 

07 = ideal (x*y -l,x +y +x+y+l,y +y +x+y+l) 

07 : Ideal of S 

18 : Id = ideal (a*x+b*y, a"2-2, b~2-3) ; 

08 : Ideal of S 

19 : ideal selectlnSubringd, gens gb Id) 

2 3 2 

09 = ideal (x - -*y ) 

2 

o9 : Ideal of S 

ilO : Ie = ideal (a*x+b*y-l, a~2-2, b~2-3) ; 
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olO 



Ideal of S 



ill 



ideal selectlnSubringd , gens gb Ie) 



oil 



4 2294 232 1 

ideal (x - 3x y + -*y - x - -*y + -) 
4 2 4 



oil : Ideal of S 

Multiplicity. The multiplicity of a zero-dimensional scheme X at a 
point p e X is defined to be the length of the local ring Ox, p - Unfortu- 
nately, we cannot work directly in the local ring in Macaulay2. What 
we can do, however, is to compute the multiplicity by computing the 
degree of the component of X supported at p\ see page 66 in Eisenbud 
and Harris (1999). 

Question 5. What is the multiplicity of the origin (0, 0, 0) as a zero 
of the polynomial equations 

x 5 + y 3 + z 3 = x 3 + y 5 + z 3 = x 3 + y 3 + z 5 = ? 

Solution. If / is the ideal generated by x 5 + y 3 + z 3 , x 3 + y 5 + z 3 and 
x 3 + y 3 + z 5 in Q[x, y, z], then the multiplicity of the origin is 



It follows that the multiplicity is the vector space dimension of the ring 



Q[x,y,z]/(p 1 (I-Q[ 

x i Vi z]{ x ,y,z)) where y?: Q[x,y,z] — > Q[x, y, z]( Xj y iZ ) is 



the natural map. Moreover, we can express this using ideal quotients: 



Carrying out this calculation in MACAULAY2, we obtain: 

11 : S = QQ[x, y, z] ; 

12 : I = ideal (x~5+y~3+z~3, x"3+y"5+z"3, x"3+y"3+z"5) ; 

02 : Ideal of S 

13 : multiplicity = degreed : saturate (I)) 

03 = 27 

Flat Families. No n- reduced schemes arise naturally as the flat limit 
of a family of reduced schemes. Exercise 111-68 in Eisenbud and Harris 
(1999) illustrates how a family of skew lines in P 3 gives a double line 
with an embedded point. 

Question 6. Let L and M be the lines in PjL, given by x = y = and 
x — tz = y + t 2 w = respectively. Show that the flat limit as t — > of 
the union LU M is the double line x 1 = y = with an embedded point 
of degree 1 located at the point (0:0:0:1). 

Solution. We find the flat limit by saturating the intersection ideal: 



dim, 



Q[x,y,z](x,y,z) 



■Q 



I ■ Q[x,y,z](x,y,z) ' 



V -\I-Q[x,y,z] {x ^ z) ) = (/:(/: {x,y,z)°°)) . 
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i2 


ci = nn r+- v ^ -7 t.t~i • 


i3 


phi = map(PP3, S, I vars PP3 ); 


o3 


RingMap PP3 < S 


i4 


L = ideal Cx, yj ; 


o4 


Ideal of S 


i5 


M = ideal (x-t*z, y-t~2*w) ; 


o5 


Ideal of S 


i6 


X = intersect (L, M) ; 


06 


Ideal of S 


i7 


Xzero = trim phi substitute (saturate (X, t) , {t => 0» 
2 2 


o7 = 


= ideal (y*z, y , x*y, x ) 


o7 


Ideal of PP3 



This is the union of a double line and an embedded point of degree 1. 

18 : use PP3; 

19 : intersect (ideal (x"2, y) , ideal(x, y~2, z)) 

2 2 
o9 = ideal (y*z, y , x*y, x ) 

09 : Ideal of PP3 

ilO : degree( ideal(x~2, y) / ideal(x, y~2, z)) 

010 = 1 

Bezout's Theorem. Bezout's theorem (Theorem 111-78 in Eisenbud 
and Harris, 1999) fails without the Cohen-Macaulay hypothesis. Fol- 
lowing Exercise 111-81 in Eisenbud and Harris (1999), we illustrate this 
in MACAULAY2. 

Question 7. Find irreducible closed subvarieties X and Y in P 4 with 

codim(X n7) = codim(X) + codim(y) and 

deg(Xnr) > deg(X) • deg(F) . 

Solution. We show that the assertion holds when X is the cone over the 
nonsingular rational quartic curve in P 3 and Y is a two-plane passing 
through the vertex of the cone. The computation is done as follows: 

11 : S = QQ [a, b, c, d, e] ; 

12 : quarticCone = trim minors (2, 

matrix-f-Ca, b~2, b*d, c}, {b, a*c, c"2, d}>) 

3 2 2 2 3 2 

o2 = ideal (b*c - a*d, c - b*d , a*c - b d, b - a c) 

02 : Ideal of S 

13 : plane = ideal(a, d) ; 

03 : Ideal of S 

14 : codim quarticCone + codim plane == codim (quarticCone + plane) 

04 = true 
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15 : (degree quarticCone) * (degree plane) 

05 = 4 

16 : degree (quarticCone + plane) 

06 = 5 

Constructing Blow-ups. The blow-up of a scheme X along a sub- 
scheme Y can be constructed from the Rees algebra associated to the 
ideal sheaf of Y in X; see Theorem IV-22 in Eisenbud and Harris 
(1999). Grobner basis techniques allow one to express the Rees al- 
gebra in terms of generators and relations. We demonstrate this by 
solving Exercise IV-43 in Eisenbud and Harris (1999). 

Question 8. Find the blow-up of the affine plane A 2 = Spec (k[x, 7/]) 
along the subscheme defined by (x 3 , xy , y 2 ) . 

Solution. We first provide a general function which given an ideal and 
a list of variables returns the ideal of relations for the Rees algebra. 

11 : blowUpIdeal = methodO ; 

12 : blowUpIdeal (Ideal, List) := (I, L) -> ( 

r := numgens I; 
S : = ring I ; 

kk := coef f icientRing S; 
n := numgens S; 
y := symbol y; 

St := kk[t, gens S , y_l..y_r, MonomialOrder => Eliminate 1]; 
phi := map (St, S, submatrix (vars St, {l..n})); 
F : = phi gens I ; 
local J; 

J = ideal apply(l..r, j -> y_j - t*F_(0, (j-1))); 

J = ideal selectlnSubringd , gens gb J); 

if (#L < r) then error "not enough variables"; 

R : = kk [gens S , L] ; 

theta := map(R, St, I vars R) ; 

theta J 

); 

Applying the function to our specific case yields: 

13 : S = QQ[x, y] ; 

14 : I = ideal(x~3, x*y, y~2) ; 

04 : Ideal of S 

15 : blowUpIdeal (I, {A, B, C}) 

2 2 3 2 

05 = ideal (y*B - x*C, x*B - A*C, x B - y*A, x C - y A) 
o5 : Ideal of QQ [x, y, A, B, C] 

A Classic Blow-up. We consider the blow-up of the projective plane 
P 2 at a point. Many related examples appear in section IV.2.2 of 
Eisenbud and Harris (1999). 

Question 9. Show that the following varieties are isomorphic. 



M2 AND SCHEMES 



9 



(a) the image of the rational map from P 2 to P 4 given by 



(b) the blow-up of the plane P 2 at the point (0:0:1); 

(c) the determinantal variety defined by the 2x2 minors of the matrix 



This surface is called the cubic scroll in P 4 . 

Solution. We find the ideal in part (a) by elimination theory. 

11 : PP4 = QQ[a, b, c, d, e] ; 

12 : S = QQ[r, s, t, A..E, MonomialOrder => Eliminate 3 ]; 

13 : I = ideal(A - r~2, B - s~2, C - r*s, D - r*t , E - s*t) ; 

03 : Ideal of S 

14 : phi = map(PP4, S, (matrix{{0, 0, 0»**PP4) I vars PP4) 

04 = map(PP4,S,{0, 0, 0, a, b, c, d, e}) 
o4 : RingMap PP4 < S 

15 : surfaceA = phi ideal selectlnSubringd, gens gb I) 



o5 = ideal (c*d - a*e, b*d - c*e, a*b - c ) 
o5 : Ideal of PP4 

We determine the surface in part (b) by constructing the blow-up of 
P 2 in P 2 x P 1 and then projecting its Segre embedding from P 5 into P 4 . 
Notice that its image under the Segre map lies on a hyperplane in P 5 . 

16 : R = QQ[t, x, y, z, u, v, MonomialOrder => Eliminate 1]; 

17 : blowUpIdeal = ideal selectlnSubringd , gens gb ideal (u-t*x, v-t*y)) 
o7 = ideal (y*u - x*v) 

o7 : Ideal of R 

18 : PP2xPPl = QQ[x, y, z, u, v] ; 

19 : psi = map(PP2xPPl, R, I vars PP2xPPl) ; 
o9 : RingMap PP2xPPl < R 

110 : blowUp = PP2xPPl / psi(blowUpIdeal) ; 

111 : PP5 = QQ[A, B, C, D, E, F] ; 

112 : segre = map(blowUp, PP5, matrix{{x*u, y*u, z*u, x*v, y*v, z*v») ; 
ol2 : RingMap blowUp <— PP5 

113 : ker segre 



ol3 = ideal (B - D, C*E - D*F, D - A*E, C*D - A*F) 

013 : Ideal of PP5 

il4 : theta = map( PP4, PP5, matrix{{a, c, d, c, b, e}}) 

014 = map(PP4,PP5,{a, c, d, c, b, e}) 
ol4 : RingMap PP4 < PP5 




a C d 
b d e 




2 



2 
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115 : surfaceB = trim theta ker segre 

2 

ol5 = ideal (c*d - a*e, b*d - c*e, a*b - c ) 

015 : Ideal of PP4 

Finally, we compute the surface in part (c) and apply a permutation 
of the variables to obtain the desired isomorphisms 

116 : determinantal = minors(2, matrix{{a, c, d},{b, d, e}}) 

2 

016 = ideal (- b*c + a*d, - b*d + a*e, - d + c*e) 

016 : Ideal of PP4 

117 : sigma = map(PP4, PP4, matrix{{d, e, a, c, b}}) ; 

017 : RingMap PP4 < PP4 

118 : surfaceC = sigma determinantal 

2 

018 = ideal (c*d - a*e, b*d - c*e, a*b - c ) 

018 : Ideal of PP4 

119 : surfaceA == surfaceB 

019 = true 

120 : surfaceB == surfaceC 

020 = true 

Fano Schemes. Our last example concerns the family of Fano schemes 
associated to a flat family of quadrics. We solve Exercise IV-69 in 
Eisenbud and Harris (1999). 

Question 10. Consider the one-parameter family of quadrics 

Q = V{tx 2 + ty 2 + tz 2 + w 2 ) C F 3 k[t] = Proj [k[t] [x, y, z,w\) ■ 

As the fiber Q t tends to the double plane Q , what is the flat limit of 
the Fano scheme Fi(Q t ) of lines lying on these quadric surfaces ? 

Solution. We first make the homogeneous coordinate ring of the ambi- 
ent projective 3-space and the ideal of our family of quadrics. 

11 : PP3overBase = QQ[t, x, y, z, w] ; 

12 : Qt = ideal(t*x~2+t*y~2+t*z~2+w~2) ; 
o2 : Ideal of PP3overBase 

We construct an indeterminate line in PL t i by adding parameters u, v 

and two points (A: B: C: D) and (E: F: G: H). The map <fi sends the 
variables to the coordinates of the general point on this line. 

13 : S = QQ[t, u, v, A. .H] ; 

14 : phi = map(S, PP3overBase, matrix{{t» I 

u*matrix{{A, B, C, D}}+v*matrix{{E, F, G, H}}) ; 

o4 : RingMap S < PP3overBase 

The indeterminate line is contained in our family of quadrics if and 
only if 4>{tx 2 + ty 2 + tz 2 + w 2 ) vanishes identically in u and v. Thus, 
we extract the coefficients of u and v. 
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15 : imageOfQt = phi Qt; 

05 : Ideal of S 

16 : coeff = (coef f icients ({1 ,2>, gens imageOf Qt) ) _1 ; 

1 3 

06 : Matrix S < S 

We no longer need the variables u and v. 

17 : Sprime = QQ[t, A. .H] ; 

18 : coeff = substitute (coef f , Sprime); 

1 3 
08 : Matrix Sprime < Sprime 

19 : Sbar = Sprime / (ideal coeff) ; 

To move to the Grassmannian over Q[t], we introduce a polynomial ring 
in 6 new variables corresponding to the minors of the matrix [ ^ f g ] ■ 
The map ip sends the new variables a, ... f to the appropriate minor, 
regarded as elements of Sbar. 

110 : PP5overBase = QQ [t , a..f]; 

111 : psi = matrix-[{Sbar_"t"}-} I substitute( 

exteriorPower(2, matrix{{A, B, C, D>, {E, F, G, H}» , Sbar); 

1 7 

011 : Matrix Sbar < Sbar 

112 : fanoOfQt = trim ker map(Sbar, PP5overBase, psi); 

012 : Ideal of PP5overBase 

We next determine the limit as t tends to 0. 

113 : fanoOfQO = trim substitute (saturate ( 

fanoOfQt, t) , {t => 0}) 

2 2 2 

013 = ideal (e*f, d*f , d*e, a*e + b*f , d , f , e , c*d - b*e + a*f, 

2 2 2 

b*d + c*e, a*d - c*f , a + b + c ) 

013 : Ideal of PP5overBase 

114 : degree( ideal(d, e, f, a~2+b~2+c~2) / fanoOfQO) 

014 = 2 

We see that Fi(Q ) is supported on the plane conic (d, e, /, a 2 + b 2 + c 2 ) 
and Fi(Q Q ) is not reduced — it has multiplicity two. 

From section IV. 3. 2 in Eisenbud and Harris (1999), we know that 
Fi(Qi) is the union of two conies lying in complementary planes. We 
verify this as follows: 

115 : fanoOfQl = trim substitute (saturate ( 

fanoOfQt, t) , {t => 1}) 

2 2 2 

015 = ideal (a*e + b*f , d + e + f , c*d - b*e + a*f , b*d + c*e, 

2 2 2 2 2 

a*d - c*f, c + e + f , b*c + d*e, a*c - d*f, b - e , 

2 2 
a*b + e*f, a - f ) 



12 



SMITH AND STURMFELS 



015 : Ideal of PP5overBase 

il6 : fanoOfQl == intersect(ideal(c-d, b+e, a-f , d~2+e~2+f ~2) , 
ideal(c+d, b-e, a+f , d~2+e~2+f~2) ) 

016 = true 

Thus, Fi(Q ) is the double conic obtained when the two conies in 
Fi{Qi) move together. 
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